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INTRODUCTION 


In recent years the world's food sltuatir/n has emphasized the need 
for timely information on world-wide crop production. Relatively few 
countries, however, have reliable methods for gathering and reporting 
crop production information. 

Remote sensing from aerospace platforms can provide information 
about crops and soils which could be useful to crop production 
forecasting systems. The feasibility of utilizing multispeotral data 
from satellites to identify and measure crop area has been demonstrated 
(1^) relatively little research has been conducted on .developing the 
capability of using multispeotral data to provide Information about crop 
condition and yield. If this spectrally derived information can be 
combined effectively with crop models which depict limitations imposed 
on crop yields by weather and climate, then potentially much better 
information about crop yield and production can be gained. 

Solar radiation is the source of energy for photosynthesis, the 
initial process that green plants use to convert carbon dioxide and 
water into simple sugars. Other plant processes convert these initial 
products of photosynthesis into dry matter including carbohydrates, 
proteins, and oils. Solar radiation is available as an energy source 

for- plants only when it interacts with leaves. In a healthy crop 

adequately supplied with water, the production of dry matter is 

proportional to the solar radiation intercepted by the canopy. Thus, 

important components of growth and yield are the amount and duration of 
plant surface available for photosynthesis 

In theory, the production of dry matter (DM) over time period t, 
beginning at emergence and ending at maturity, can be related to the 
proportion (P) of the incident light (SR) intercepted by the crop using 
the following equation from Steven (1902): 

m 

DM = J E P SF dt. [1] 

e 

E is the efficiency of conversion of solar energy into dry matter and 
typically ranges from 1 to 3 g/M Joule (15). This equation can be used 
to predict dry matter production if P is known. Current methods to 
measure the radiation intercepted by crops are laborious and limit use 
of such crop models to small research plots. If the proportion of 
energy available for crop growth could be estimated reliably using 
multispeotral satellite data, then the capability to estimate crop 
production for large regions should be improved significantly. 

In practice, although solar radiation is essential for 

photosynthesis, it is only one of several factors interacting to 
Influence crop yields. Other factors essential to crop growth and yield 
are water, temperature, nutrients, and carbon dioxide. Any serious and 
comprehensive effort to estimate crop yields also must assess the impact 
of these other factors. 
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Our overall objeobive is to develop methods for combining spectral 
and meteorological data in crop yield models which are capable of 
providing accurate estimates of crop condition and yields throughout the 
growing season, This paper presents cur initial bests of these concepts 
using spectral and agronomic data acquired in controlled experiments. 
Future research will extend those methods that best estimate yields at 
agricultural experiment stations bo large areas using spectral data 
acquired by satellites. 


MATERIALS AND METHODS 
Design of Experiment 

Spectral and agronomic data used in these analyses were acquired at 
the Purdue University Agronomy Farm in 1979 and 198O. A full season 
corn ( Zea mays L.) hybrid, Becks 65X, was grown on two soil types. When 
dry, the Chalmers silt loam (Typic Argiaquoll) had a dark gray (lOYR 
VI) surface while the Toronto silt loam (Odollic Oohraqualf) had a 
light gray (lOYR 6/1) surface. These two soils were spectrally distinct 
in both visible and near infrared reflectance factor (11). 

Prior to planting 200, 50, and 95 Rg/na of N, P, and K respectively 
were applied uniformly to both soils bo minimize the risk that corn 
growth might be limited by nutrient availability. Daily meteorological 
data were recorded at the cooperative National Weather Service station 
(West Lafayette 6NW) which was within 200 ra of the fields. Incoming 

solar radiation was measured with an Eppley Precision _ Spectral 

Pyranoraeter and recorded as total Langleys/day ( cal cm""^ day Daily 
maximum and minimum air temperatures were measured with liquid -in -glass 
thermometers in a standard shelter. Soil moisture in the top 105 cm 
(depth of drainage tiles) was estimated using a soil moisture balance 
model (19). 

Within each soil type two completely randomized blocks with three 
plant populations (25,000, 50,000, and 75,000 plants/ha) were planted in 
76 cm wide north-south rows on 2, 16 and 30 May 1979 and 7, 22 May and 

11 June 1980. Additional plots of the 50,000 plants/ha treatments were 
planted on 16, 29 May, I8 June, and 3 July 1980, These treatments 
represented a wide range of planting dates and plant populations 
expected in corn fields in Indiana. 


Canopy Characterization 

Agronomic variables measured and observed coincidentally with 
reflectance factor measurements included; plant height, leaf area index 
(LAI), development stage (8), total fresh and dry biomass, dry stalk 
(including leaf sheath), ear, and green leaf blade (lamina) weights. 
Percent soil cover (defined as the percentage of soil covered by 
vegetation) was determined by placing a grid over a vertical photograph 
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and counting the interscotiona occupied by green vegetation. Crop 
biomass was estimated by harvesting three plants in 1979 and four plants 
in 1980 from each plot. Each sample was weighed immediately (fresh 
weight), separated into its components, dried at 75® C and reweighed. 
The area of a random subsample of green leaf blades from each plot was 
measured with an electronic area meter (LI-COR, Model LI-3000) and the 
green leaf area to leaf dry weight ratio was calculated. Leaf area 
index was calculated using this leaf area/weight ratio, the total dry 
weight of green leaves from all plants sampled, and the soil area 
represented. Grain was harvested by hand from the center four rows of 
each plot (11.6 m^) , dried, weighed, and corrected to 15.5 percent 
moisture. Visual assessment of the soil moisture and crop oondition 
were made during the spectral data collection. Crop oondition 

assessment included evaluations of lodging and hail and insect damage. 


Spectral Measurements 

Radiance measurements, used to determine reflectance factor, were 
acquired with a Landsat-band radiometer (Exotech 100) throughout the 
growing season in each year. Robinson and Biehl (1979) describe the 
conditions and procedures for obtaining the RF data. The Exotech 100 
has a 15 degree field of view and acquired data in the following 
wavelength regions: 0.5-0. 6, 0.6-0.7i 0,7-0. 8, and 0. 8-1.1 ym. Data 

were taken only when there were no clouds over or in the vicinity of the 
sun and vdien the solar elevation was at least 45 degrees above the 
horizon. 

The radiometer and a 35 -mm camera were attached to a boom mounted 
on a piokup truck and elevated 5.2 m above the soil in 1979 and 7.6 m in 
1980. After the instruments were leveled for a nadir look angle, two 
measurements were taken - one centered over the row and one centered 
between rows - in each plot to better estimate the overall canopy 
response (5). A color photograph which included the ar-ea viewed by the 
radiometer was taken vertically over each plot and used to determine 
percent soil cover. 


Analysis 

Two methods of estimating the proportion of solar radiation 
intercepted by corn canopies were examined. First, the proportion of 
intercepted radiation (SRIm) was described as a function of measured LAI 
using the following equation from Linvill et al. (1978): 

SRI^ = [1 - exp(-0.79 LAI)]. [2] 

This is an application of Bouguer's Law using LAI and an extinction 
coefficient of -0.79. When LAI is 0, no energy is Intercepted. When 

LAI is 2.8, approximately 90 % of the visible solar radiation is 
intercepted by the canopy and is potentially useful to the crop. 
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The aeoond tnebhod esbimafcee SRI as a fUnobion of the spectral 
variable greenness using the following equation; 

SRig s -0.1613 + 0.0811 G - 0.0015 G^- [3] 

where G Is the green vegetation index or ’’greenness” for refleotanoe 
factor data (13). Greenness was calculated as follows; Greenness » 
(-0.H89^ RPj+ -0,6125 RF?^ + 0,1739 RFj 4- 0.5953 RF4), where RR to RF,, 
refer to the reflectance factor in each of four bands of the radiometer. 
This spectrally estimated proportion of radiation intercepted is called 
SRIg bo distinguish it from which is estimated using measured LAI, 

Factors involved in determining the interception of radiation by 
vegetation include solar azimuth and zenith angles, spectral properties 
of canopy elements, leaf area index, leaf angle distribution, leaf size 
and shape, and leaf movement due to wind, wilting, jmid phototropism 
( 16 ), The concept of a simple exponential extinction of radiation 
appears to be widely applicable. Norman (1980) presents a summary of 
extinction coefficients reported for various canopies and sun elevation 
angles. For corn the extinction coefficients ranged from -1.5 for low 
sun angles (10“) to -O .56 for high solar elevation angles (70°). The 
extinction coefficient used in this research is within the range of 
values presented by Norman for approximately iJ5 degree solar elevation. 
Additional research is needed to oharaoberize and model the changes in 
extinction coefficient if this approach is to be used quantitabively, 

SRIg predicted as a function of greenness and SRlj,j predicted as a 
function of measured LAI were calculated for each day that appropriate 
spectral and agronomic data were acquired and linearly interpolated for 
intermediate days throughout the growing season for each plot. Daily 
values of SRIg and SRI^, were accumulated from planting to physiological 
maturity. Direct correlations of final grain yields with these 
accumulated indexes were examined. 

To account for variability in temperature and plant water status, 
the performances of spectrally estimated SRIg and measured SRIm 
were compared using the Energy-Crop Growth (ECG) model (3,^) which 
combines the concept of intercepted solar radiation with a moisture 
stress term and a temperature function. The ECG model used was; 

mature 

ECG„j s Z (SRi/600)(SRl„j )(WF^)(FT^) [4] 

implant ^ 

where SR is the daily solar radiation in__gal cm ' day and 600 is the 
approximate latent heat of water (cal om""^), HF is the ratio of daily 
evapo transpiration to potential evapotranspiration (19), and FT is a 
daily tenperature function (3). For the spectral ECGg model SRIg was 
substituted directly for SRI^ in equation [i|J. 
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RESULTS AND DISCUSSION 

Relation of Canopy Refleobanoe to LAI and SRI 

Leaf area index (LAI) and the proportion of solar radiation 
intercepted (SRIg) by these corn canopies were described as functions of 
several spectral variables and transformations using regression 
analyses. Previous research has indicated that greenness is highly 
correlated to LAI and percent soil cover but relatively insensitive to 
soil color (11,20) in each year and for the combined data greenness 
predicted SRI better (higher R^) than LAI, 

The relationships of LAI and SRI to greenness for all plots of 
corn in 1979 and I98O are shown in Figure 1. Undoubtedly planting date, 
plant population, and soil type contt.''ibuted to the scatter about the 
regression line, but when included aii terms in the regression models, 
they contributed very little additional information. Errors in 
measuring LAI also account for a portion of the scatter. Nevertheless, 
LAI and SRIg predicted as a function of greenness permit the concepts 
and models developed at the Purdue Agronomy Farm to be used in 
situations where only spectral data may be available. 

The response of greenness to LAI appears asymptotic for LAIs 
greater than approximately 5.0. This is consistent with infinite LAI 
refleotanoe of single leaves where visible reflectance vias minimized 
with two layers of leaves and near infrared reflectance was maximized 
with six to eight layers (6). Furthermore, since most of the incoming 
radiation is intercepted by canopies with LAIs less than 4.0 the 
importance of accurately estimating LAIs greater than 4.0 is greatly 
diminished. The critical value of LAI - when 90 percent or more of the 
radiation is intercepted - increases as the extinction coefficient 
decreases. Factors affecting the extinction coefficients of crop 

canopies have been described (2, 16). 


Relation of Spectral Variables to Yield 

Corn is a determinant crop; that is, corn completes its vegetative 
development, producing all of its leaves, then shifts to reproductive 
development, producing and filling its grain. If corn grain yields are 
strongly related to the total amount of green leaf area present, then 
knowing the maximum LAI or maximum greenness of corn should be an 
important predictor of grain yields. Figure 2 illustrates the errors of 
such logic. Grain yields are plotted as a function of maximum values of 
greenness which occurred near tasseling. The poor relationship shown in 
Figure 2 emphasizes the limited value of a single observation of 
greenness for predicting corn yields. 

A more important and useful indicator of grain yields is the 
seasonal duration of LAI and not simply the maximum LAX achieved. 
Previous researchers have reached similar conclusions for other crops 
(4,12,18). Two variables which represent the integral of LAI over time 
during the growing season are accumulations of the daily values of SRIj^ 
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Figure 1. Relationships of leaf area index (LAI) solar radiation 
intercepted (SRIg) to the spectral variable greenness. LAI = 
(-0.1364 + 0.064? Greenness + 0.004? Greenness^). SRIg = (-0.1613 + 
0.0811 Greenness - 0.00 15 Greenness^). 
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Figure 2. Corn, grain yields as a f\inotion of maximum greenness. 




and SRIg. These summed values of SRIm and SRIg represent the 
proportions of the solar radiation impinging on the corn field during 
the growing season that were intercepted by the oorn canopy and thus 
were potentially available for photosynthesis (J?’'igure 3)* Both 
estimates of intercepted solar radiation were assooiated with 
approximately 65 percent of the variation in grain yields. 

One problem in crop response to solar radiation is the confounding 
of solar radiation, temperature, and plant moisture stress on plant 
growth and yields. Thus applicability of the theoretical model 
presented in Eq. fl] is extremely limited. Several researchers have 
demonstrated for various crops that the reduction in crop growth is 
proportional to the reduction in evapotranspiration (ET) from potential 
evapotranspiration (PET) (4,7»9,10). Coelho and Dale ( 198 O) combined 
temperature and moisture stress with an estimate of intercepted solar 
radiation using measured LAI in an Energy-Crop Growth (ECO) model, Eq. 
[4). They used this model to simulate the daily effects of weather on 
oorn growth. 

The sum of the daily values of ECO using SRI or SRI to estimate 
intercepted solar radiation are plotted with grain yields in Figure 4. 
Both of the ECG models have slightly larger coefficients of 
dat'd rmlnation (r^) than the SAIs and SRlm models. Although the ten 
planting date-years represented a range of temperature and day length 
regimes, significant water stress was not evident among the treatraents>> 
In this ease, both of the intercepted radiation variables (SRIh and SRIn) 
alone were nearly as highly correlated to yields as the meteorological 
(ECG) models. In situations where moisture or temperature limits yield, 
it is postulated that either of the two ECG models should be superior to 
the intercepted solar radiation models for predicting oorn yields. 

The concept of combining spectrally derived estimates of SRIs with 
meteorological models should permit implementation of crop growth and 
yield models for large areas. A spectral-meteorological system of crop 
forecasting could exploit the frequent temporal sampling of weather data 
(e.g. dail,v or hourly) with the high spatial resolution typical of earth 
observing satellites (e.g. Landsat MSS). 

An example of this concept showing how multispectral data from 
satellites could be used to estimate crop yields is presented in Figure 
5 . Percent soil cover and leaf area index strongly influence the 
reflectance of radiation from crop canopies (11,20). Estimates of 
percent soil cover and leaf area index, obtained from Landsat MSS data 
along with daily solar radiation from ground stations or meteorological 
satellites (1) will permit calculation of the amount of radiation 
intercepted by crops. This intercepted solar radiation variable 
combined with observed temperature and precipitation data and integrated 
over the growing season should account £or much of the variation in corn 
yields. The multispectral data from satellites would form the basis for 
estimating crop growth and yields over regions where ground observations 
may be difficult or impossible to obtain. 
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ITlgure 3. Corn grain yields as a function of the accumulated 
proportions of solar radiation intercepted by the corn canopies. 
SRIj„ is estimated using measured LAI and SRIg is estimated using 
spectral data. The symbols are relative planting dates for the 
two years with "0" representing the first planting date in 1979 
and "9" representing the last planting date in 1980, 
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SUMMARY AND CONCLUSIONS 


Spectral and agronomic data for corn were acquired throughout two 
growing seasons. Two variables, SRI,„ and SRIg , representing the 
proportions of solar radiation potentially available for photosynthesis 
were estimated. SRI„^ was calculated by Bouguer’s Law using measured LAI 
while SRIg was estimated using the spectral variable of greenness. 

Various methods or models for predicting yields were examined in 
this research. One method using a single observation of the maximum ^ 

greenness had the lowest correlation with corn grain yields. Methods 
which integrated spectral or meteorological variables over the growing 
season (e.g,, SRIm, SRIg) accounted for over half of the variation in 
grain yields. Models which simulated the daily effects of weather on 
corn growth (e.g., ECGn^, ECGg) had the highest correlations to grain 
yields. 

We conclude from this research that the concept of estimating 
intercepted solar radiation using spectral data represents a viable 
approach for merging spectral and meteorological data in crop yield 
models. This scheme is consistent with the spectral-physiological 

modeling approaches proposed by Wiegand et al. (1979) but may avoid some 
of the problems associated estimating LAI directly from spectral data. 

This concept may be extended to large areas using Landsat MSS data to 
estimate SRIg for as many fields as are of interest. We are currently 
assembling the necessary data to evaluate this concept using Landsat MSS 
data from commercial fields in the U.S. Corn Belt. 
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